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MFTHOD FOR DETER MINING SUSCEPTIBILI TY TO BONE DAMAGF BY SCREENING POLYMORPHISMSJN_THE 

VTTAMIN D RECEPTOR GENE 

Ttwpresent invention relates to a prognostic method and prognostic means based on 
polymorphisms in the vitamin D receptor and collagen loci genes. In particular, the 
5 present invention relates to a method for determining susceptibility to bone damage by 
screening for polymorphisms in vitamin D receptor or collagen led genes. 

Osteoporosis is a common disease characterized by reduced bone mineral density 
(BMD), deterioration of bone micro-architecture and increased risk of bone damage, 

10 such as fracture.' It is a major public health problem which affects quality of life and 
increases costs to health care providers. In European populations, one in three women 
and one in twelve men over the age of fifty is at risk. The disease effects 25 million 
people in the USA, where the incidence of disease is 25% higher than it is in the UK, 
and a further 50 million people in Japan and Europe combined. It is estimated that by 

15 the middle of the next century the number of osteoporosis sufferers will double in the 
West, but may increase six-fold in Asia and South-America. Fracture is the most 
serious endpoint of osteoporosis, particularly fracture of the hip which affects up to 
1.7 million people worldwide each year. It is estimated that by the year 2050, the 
number of hip fractures worldwide will increase to over 6 million, as life expectancy 

20 and age of the population increase. 

Treatment of osteoporosis is unsatisfactory. In particular, once bone damage has 
occurred as a result of osteoporosis, there is Utile a physician can do other than let the 
bone heal. In the elderly, this may be a slow and painful process. Diagnosis of those 
25 at risk of developing osteoporosis allows more effective preventative measures. 
Strategies for the prevention of this disease include development of bone density in 
early adulthood, and minimisation of bone loss in later life. Changes in lifestyle, 
nutrition and hormonal factors have been shown to affect bone loss. 7 " 14 
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Osteoporosis can be considered a complex genetic trait with variants of several genes 
underlying the genetic determination of the variability of the phenotype. Low bone 
mineral density (BMD) is an important risk factor for fractures, the clinically most 
relevant feature of osteoporosis. Segregation analysis in families has shown that BMD 
is under polygenic control 12 ' 13 while, in addition, biochemical markers of bone 
turnover have also been shown to have strong genetic components 14 " 16 . Several 
candidate genes have been analysed in relation to BMD but the most widely studied 
gene in this respect, the vitamin D receptor (VDR) gene, explains only a small part of 
the genetic effect on BMD S . Numerous studies, focussing on the Bsml allele of the 
vitamin D receptor gene have concluded that absence of the restriction site correlates 
with low bone mineral density. Most genetic analyses have focussed on BMD as a 
determinant of fracture risk and not so much on fractures themselves as an endpoint in 
the analysis. Recently, an Spl polymorphism in the COLIAl gene encoding the most 
abundant bone matrix protein, was found to be associated with reduced BMD and, 
more importantly, also with increased risk of osteoporotic fracture 9 ' 10 . An emerging 
theme from these studies seems to be the association of the absence of the Bsml 
restriction site with reduced bone mineral density, which in turn signifies increased 
risk of bone damage such as fracture. 

Few studies have addressed genetic association with the clinically most important 
endpoint of osteoporosis, namely bone damage, in particular fracture. Accordingly, it 
is an object of the present invention to improve the prognosis of predisposition or 
susceptibility to bone damage. 

In a first aspect of the present invention, there is provided a method of determining 
susceptibility to bone damage in a subject, said method comprising analysing genetic 
material of a subject to determine which allele(s) of the vitamin D receptor gene is/are 
present. 
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Thus, the present invention satisfies the pressing need for identification of those 
individuals susceptible to bone damage, thus facilitating the development of 
preventative measures. For example, those at risk may avoid damage by modifying 
their lifestyle and implementing bone strengthening measures, such as regular 

5 exercise and a healthy diet. Typically, the method of the present invention 
comprises analysis of polymorphisms in the VDR gene to determine susceptibility to 
bone damage. The method may include determining whether one or more particular 
alleles are present, or which combination of those alleles (i.e. the haplotype) is 
present. The method may further comprise determining whether subjects are 

10 homozygous or heterozygous for alleles or haplotypes of the vitamin D receptor 
gene. 

Vitamin D is a potent regulator of bone and calcium homeostasis, as well as of cellular 
differentiation and replication in many tissues, and mediates its effects .through the 

15 vitamin D receptor (VDR). Cloning of the vitamin D receptor has shown it to be a 
member of the ligand-activated superfamily, which are natural regulators of a number 
of physiological and developmental processes. Evidence suggests that the vitamin D 
receptor activates expression of the osteocalcin gene through interaction with a 
palindromic sequence in the promoter of the gene. 27 The osteocalcin gene product is a 

20 marker of bone turnover in normal and disease states, and inter-individual variation in 
its circulating levels have been associated with polymorphisms in the vitamin D 
receptor gene. 

The vitamin D receptor gene (12ql2) comprises inherited polymorphisms between 
25 exon 7 and the 3* UTR of the VDR gene, as shown in Figure 1. These alleles are 
denoted B/b, A/a and T/t for restriction enzyme sites BsmI, Apal and Taql respectively 
(or enzymatic or chemical procedures having similar specificity), where a lower case 
letter denotes the presence of a wild type restriction site which is capable of being 
cleaved, and a capital letter denotes the presence of a mutant restriction enzyme site 
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which is not capable of being cleaved by the relevant restriction enzyme. For the 
purposes of the present invention, determination of which alleles are present in a 
particular gene may be referred to as determining the genotype of a subject for a 
particular gene. It is apparent from the above that each copy of the vitamin D receptor 
gene will comprise a specific combination of the three alleles, this combination being 
referred to as the haplotype of the gene. For example, the haplotype may be baT, 
indicating the presence of cleavable Bsml and Apal sites, and a non-cleavable Taql 
site. Direct haplotyping of the VDR gene has allowed five different haplotypes to be 
determined, of which three are common. 11 

The present invention is based upon the surprising observation of a correlation 
between the presence of the b allele of the vitamin D receptor and susceptibility to/or 
risk of (where the terms are used interchangeably) bone damage, such as fracture. The 
invention goes further to show that presence of the a and/or T alleles, and in particular 
the haplotype baT is/are associated with increased risk of bone damage. A subject 
having the baT haplotype may show a higher risk of fracture compared to a subject 
having the bAt, or BAt haplotype which confers the lowest risk of fracture. The 
results are unexpected, as previous studies have shown the b allele, and particularly 
the baT haplotype, to be associated with high bone density, which itself is not 
associated with fracture risk. Thus, in contrast to previous results, the present 
inventions have shown that susceptibility to bone damage is independent of bone 
mineral density in a subject. By screening for the presence of alleles of the vitamin D 
receptor gene, susceptibility to bone damage may be assessed without the need for 
analysis of bone mineral density. 

Preferably, the method of the first aspect of the present invention further comprises 
determining whether the alleles present are associated with risk of bone damage. This 
may be performed by comparing the alleles present in a subject with those known to 
be associated with risk of bone damage. For example, a visual aid detailing alleles and 
the relative risk of bone damage associated therewith may be used to determine 
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whether the genotype of the subject is associated with a high or low risk of bone 
damage. 

The first aspect of the present invention may also comprise the additional step of 
determining aspects of calcium metabolism, such as calcium levels in a subject. 
Preferably, the daily calcium intake is measured. This feature of the first aspect is 
based on the observation that the correlation between vitamin D receptor genotype 
and bone damage may be dependent upon dietary calcium intake. Specifically, in 
subjects having low calcium intake, genotype dependent risks may be greater. 

The method of the present invention may be performed in vitro. Preferably, the 
method is performed on tissue or fluid removed from the body of the subject. Thus, 
the present invention relates to a non-invasive diagnostic method, the results of which 
provide an indication of susceptibility to bone damage but do not lead to a diagnosis 
upon which an immediate medical decision regarding treatment has to be made. 

The present invention may be performed on any subject for whom it is desirable to 
determine risk of bone damage. Preferably, the subject may be a mammal. Most 
preferably, the subject is a human, preferably a female. 

Bone damage may be any form of structural damage including fractures, breaks, or 
chips. The term may also include biological degradation or deterioration of bone. 
Typically, the term bone damage does not include low bone density. This is in line 
with the finding that risk of bone damage is independent of bone density. Fracture 
may be defined as the clinically most important endpoint, and thus the method of the 
first aspect of the invention preferably relates to a method of determining risk of 
fracture. Although such bone damage will usually be the result of osteoporosis, it is 
irrelevant for the purposes of the present invention whether a subject has first been 
diagnosed as having osteoporosis. 
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In a preferred feature of the first aspect of the present invention, there is provided a 
method of determining susceptibility to bone damage, said method comprising 
analysing the genetic material of a subject to determine which of the B/b, A/a or T/t 
alleles of the vitamin D receptor are present. The method may comprise determining 

5 whether more than one of the above alleles is present. The subject may be further 
classified as heterozygous or homozygous for one or more of these alleles. Preferably, 
the method comprises the additional step of determining whether the aliele(s) present 
are associated with risk of bone damage, wherein presence of the b or a alleles is 
associated with increased risk of bone damage, and presence of the t allele is 

10 associated with reduced risk of bone damage. Homozygosity for the a or b allele may 
further increase the susceptibility to bone damage in a subject, while homozygosity for 
the t allele may further decrease susceptibility. 

In a preferred feature of the first aspect of the present invention, there is provided a 
15 method of determining susceptibility to bone damage in a subject, said method 
comprising analysing the genetic material of a subject to determine the haplotype of 
the Bsml, Apal and Taql alleles at the vitamin D receptor. Preferably, said method 
comprises determining whether the haplotype of the subject is associated with risk of 
bone damage, wherein the haplotype baT is associated with high risk of bone damage. 
20 A subject homozygous for said haplotype may be at a higher risk of bone damage than 
those heterozygous for the haplotype. 

In a preferred feature of the first aspect, there is provided a method of determining 
susceptibility to bone damage, said method further comprising analysing the genetic 
25 material of a subject to determine which allele of the collagen led (17q22) gene is 
present. An allele in this gene may be denoted by S/s, where s indicates the presence 
of a G to T polymorphism at the Spl restriction site. This feature of the present 
invention is based on the observation that the correlation between increased fracture 
risk and VDR genotype may be collagen Ial genotype dependent. Specifically, in 
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those subjects having vitamin D receptor alleles indicating risk of bone damage, the 
risk may be increased in those subjects also carrying a G to T polymorphism at the Spl 
site of the collagen Ial gene. 

In a preferred feature of the first aspect, there is provided a method of determining 
susceptibility to bone damage, said method comprising determining the copy number 
of the B/b, A/a or T/t alleles or haplotype of the vitamin D receptor gene and/or the 
S/s alleles of the collagen Ial gene, where an increase in copy number is associated 
with increased risk of bone damage. 

The present invention may be performed using any suitable method known in the art. 
Preferably, a tissue or fluid sample is first removed from a subject. Examples of 
suitable samples include blood, mouth or cheek cells, and hair samples containing 
roots. Other suitable samples would be known to the person skilled in the art. The 
genetic material is then extracted from the sample, using any suitable method. The 
genetic material may be DNA or RNA, although preferably DNA is used. For 
example, the DNA may be extracted using the technique described in Sambrook et al 
(Molecular Cloning - A Laboratory Manual, Cold Spring Harbor Laboratory Press). 
Determination of the genotype of a subject may then be carried out using the extracted 
DNA, employing any one of the following techniques: 

■ Southern blot analysis following digestion with one or more appropriate restriction 
enzymes. 

• PCR amplification followed by digestion with one or more appropriate restriction 
enzymes and, optionally, separation of digestion products by gel electrophoresis. 

■ Sequencing of a relevant gene fragment by any suitable method. 

■ Visualization of heteroduplex patterns, for example on PAA or agarose gels, where 
different patterns may indicate the presence of one or more specific alleles. 

■ Separation of DNA fragments using denaturing gradient gels, wherein the degree 
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of separation will depend upon the presence or absence of one or more 
polymorphic restriction sites. 

■ Separation using SSCP analysis, the patterns of which will depend upon the 
presence or absence of one or more polymorphic restriction sites. 

■ Use of allele specific oligonucleotides, hybridization patterns of which will be 
specific for various combinations of alleles. 

■ Methods such as OLA, Taqman or dot-blot for the detection of known mutations. 

■ Visualization of DNA sites using fluorescent labelled probes for alleles of interest. 

■ RFLP analysis. 



Where it is desirable to use particular restriction enzymes in performing the present 
invention, the skilled person will understand that enzymatic or chemical procedures 
having similar specificities may also be used. For example, restriction enzymes 
having similar specificity (isoschizomers) to those described herein may be used, or 
15 chemical degredation procedures with DNA or RNA cutting specificity. 

Other techniques suitable for determining the genotype of a subject may be used in the 
present invention. 

20 Where the haplotype of a gene is to be determined, it is preferable to use a direct 
haplotyping method, as described in Uitterlinden et al n . In such a method, the 
relevant portion of the gene is amplified and then subjected to restriction enzyme 
digestion, in order to determine the presence or absence of restriction enzyme sites. 
Thus, for example, where the haplotype of the vitamin D receptor gene is to be 

25 determined, the portion of the gene between exon 7 and the 3 ' UTR may be amplified, 
and the amplified DNA digested with the Bsml, Apal or Taql restriction enzymes. Gel 
analysis may then be used to determine which alleles are present. 



30 



Preferably, a fragment may be amplified using polymerase chain reaction (PGR) 
techniques, to produce copies which, where the fragment is of the vitamin D receptor, 
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are at least 1000 base pairs in length, and most preferably at least 1800 base pairs in 
length. Where the fragment to be amplified is of the collagen Ial gene, PCR primers 
may be selected to amplify a fragment which is at least 50 base pairs in length, 
preferably at least 200 base pairs in length. PCR techniques are well known in the art, 
5 and it is within the ambit of the skilled person to identify primers for amplification of 
the appropriate region of the above genes, namely the region from exon 7 to the 3' 
UTR of the vitamin D receptor gene and the first intron of the collagen Ial gene. PCR 
techniques are described in EP-A-0200362 and EP-A-0201 184. 
I 

10 In a preferred feature of the first aspect, there is provided a method of determining 
susceptibility to bone damage in a subject, said method comprising amplifying a 
fragment comprising a portion of the region from exon 7 to the 3* UTR of the vitamin 
D receptor gene, and determining which allele(s) in the vitamin D receptor is/are 
present. Primers suitable for amplification of said portion of the vitamin D receptor 

15 gene would be readily available to a person skilled in the art. Examples of such 
primers include: 

S^eeAAeTeCTCKTGGeTGAGGXCIG-y 

20 

In features of the present invention where it is desirable to determine which allele of 
the collagen Ial gene is present, at least a portion of the first intron of the collagen 
Ial gene may be amplified, followed by determination of the presence of a Spl 
restriction site. Suitable primers include: 

25 

/tfr. ^^AAC^CXGGAGIATnUCGQAC 

-n) U ^ereeAeeee?^TGGT-GGee^ 



Additional primer sequences are described in Grant et af. 
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Where the amplified portion of the gene is larger than the above defined portions of 
the genes containing the relevant alleles, it is preferable to avoid the inclusion of 
vitamin D receptor or collagen Ial gene sequences which comprise any one of the 
Bsml, Apal or Taql restriction sites of the vitamin D receptor gene, or the Spl site of 
the collagen Ial gene. 

In a second aspect of the present invention, there is provided a method of therapy, said 
method comprising treating a subject diagnosed as being at risk of bone damage, to 
reduce the risk of bone damage. Preferably, the subject is diagnosed as being at risk of 
bone damage in accordance with the first aspect of the present invention. 

Therapy may in the form of preventative or palliative care. Suitable treatments 
include modifications to lifestyle, regular exercise and changes in diet to strengthen 
bones, and hormone therapy. Suitable treatments, including pharmaceutical 
preparations to reduce bone loss, would be known to physicians and persons skilled in 
the art. Examples include anabolic steroids, bisphosphonates, vitamin D preparations, 
calcium supplements and Hormone Replacement Therapy. 

In a third aspect of the present invention, there is provided a method of predicting the 
response of a subject to treatment, said method comprising analysing genetic material 
of a subject to determine which allele(s) of the vitamin D receptor gene and/or 
collagen Ial gene is/are present. Preferably, the method includes determining whether 
the subject is susceptible to bone damage. Where a subject has been determined as 
susceptible to bone damage, the method may further comprise administering the 
appropriate treatment. The effect of a therapeutic or preventative agent may depend 
on the underlying cause of the heart disease, and in some cases it may be preferable to 
avoid the use of certain treatments. This aspect of the present invention may also be 
useful for identifying agents which may be used in the treatment of bone damage. 
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In a fourth aspect of the present invention, there is provided use of a kit to determine 
which allele(s) of the vitamin D receptor gene and/or collagen Ial gene is/are present, 
said kit comprising (i) one or more nucleic acid primer molecules for amplification of 
a portion of the vitamin D receptor and/or collagen Ial genes, and (ii) means for 
determining which allele(s) is/are present in said genes. 

Preferably, the primer molecules are suitable for amplification of at least a portion of 
the region between exon 7 and the 3TJTR of the vitamin D receptor gene, and/or a 
portion of the first intron of the collagen Ial gene. Examples of suitable primers are 
described above. 

Means for determining which allele(s) is/are present in the vitamin D receptor gene, 
and/or collagen Ial gene may include any reagents or molecules necessary for use in 
any of the methods described above. For example, where PCR followed by DNA 
digestion is used, said means preferably include PCR reagents and one or more of the 
Bsml y Apal, Taql or Spl restriction enzymes. Where the method employs Southern 
Blotting, heteroduplex visualization, or fluorescent labelling techniques for example, 
probes which bind to the appropriate regions of the vitamin D receptor gene, and/or 
collagen Ial gene may be included. Where necessary, such probes may be labelled to 
allow detection, for example by nick-translation, radio- or fluorescent-labelling, or 
random primer extension whereby the non-labelled nucleotides serve as a template for 
the synthesis of labelled molecules. Other methods of labelling probes are well known 
in the art. 

In a preferred feature of the fourth aspect of the present invention, there is provided 
use of a kit further comprising means for indicating correlation between the genotype 
of a subject and risk of bone damage. Said means may be in the form of a chart or 
visual aid, which indicate that presence of the b allele or baT haplotype, and the S 
allele of the collagen Ial gene is associated with increased fracture risk. 
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In a fifth aspect of the present invention, there is provided a kit for determining risk of 
bone damage in a subject, said kit comprising (i) one or more nucleic acid primer 
molecules for amplification of a portion of the vitamin D receptor and/or collagen Iccl 
genes; (ii) means for determining which allele(s) is/are present in said genes; and (iii) 
means for indicating correlation between the allele(s) and risk of bone damage. 

In a preferred feature of the fifth aspect, the kit may also comprise DNA control 
samples, for comparison with DNA sequences of a subject. The control samples may 
comprise the sequence of one or more alleles of the vitamin D receptor or collagen Ial 
genes, or may comprise the sequence of various haplotypes. 

Preferred features of each aspect of the present invention are as for each other aspect, 
mutatis mutandis. 

The present invention will now be described in detail with reference to the following 
examples and figures in which: 

FIGURE 1 is a schematic presentation o^the region between exon 7 and the 3 * UTR of 
the vitamin D receptor gene. / 

20 

Example 1 - Interaction between the vitamin D receptor gene and collagen type 
I 1 gene in determining susceptibility for osteoporotic fracture in 
postmenopausal women 

25 In a large population-based study 97 fracture cases were recorded among 1004 
postmenopausal women aged 55-80 years. We found a VDR haplotype, constructed 
from three adjacent restriction fragment length polymorphisms (RFLPs) 11 to be over- 
represented among fracture cases (p=0.009) corresponding to a Relative Risk of 1.8 
(95% confidence interval 1.0-3.3) for heterozygous carriers and 2.6 (95%CI 1.4-5.0) 
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for homozygous carriers of the risk haplotype. The effect was similar for vertebral and 
nonvertebral fractures and, most importantly, independent of BMD. We observed 
significant interaction (p=0.03) between VDR and COLIA1 genotype effects. Fracture 
risk was not VDR genotype-dependent in the COLIA1 "reference" group (genotype 
GG) while in the COLIA1 "risk" group (genotypes GT and TT) the relative risk of 
fracture was 2.1 (95%CI 1.0-4,4) for heterozygous and 4.4 (95%CI 2.0-9.4) for 
homozygous carriers of the VDR risk haplotype. We conclude that both the VDR and 
COLIA1 polymorphisms are genetic markers for osteoporotic fracture in women, 
independent of BMD. Our data indicate that interlocus interaction is likely to be an 
important component of osteoporotic fracture risk. 

We first analysed the relation between fractures and VDR genotype and, second, 
studied interaction between the VDR and the COLIA1 polymorphisms. 

For COLIA1 we investigated the Spl G to T polymorphism 9 in relation to haplotypes 
constructed of three adjacent restriction fragment length polymorphisms (RFLPs) in 
the VDR gene, i.e. for BsmI, Apal, and TaqI, respectively 11 . When we analysed the 
distribution of fractures in women grouped according to VDR genotype we observed 
an overrepresentation of fractures in women carrying the haplotype 1 (Table 1). 
Women were subsequently grouped according to carrier status for this VDR haplotype 
as heterozygous carriers (including the genotypes 12 and 13) and homozygous carriers 
(consisting of genotype 11) of the risk haplotype and compared to women not 
carrying the haplotype (including genotypes 22, 23, and 33). No significant differences 
in known risk factors for osteoporosis could be observed between women grouped 
according to VDR haplotype 1 (Table 2). Similar results were obtained when the 
women were grouped according to VDR haplotypes 2 or 3 (data not shown). 

We then went on to determine! the distribution of fractures in women according to their 
carrier status for VDR hAplotype 1 (Table 3a). Significantly more women 
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heterozygous for VDR haplotype 1 had fractures than the women in the reference 
group and for women homozygous for the VDR haplotype 1 this difference further 
increased. When women were grouped according to VDR haplotype 2, we observed an 
under-representation in I fracture cases (p=0.002) while for VDR haplotype 3 no 
differences were observed (p=0.65;data not shown). Logistic regression analysis 
showed that women heterozygous for the VDR haplotype 1 had 1.8 times the risk for 
fractures compared to women in the reference group. This was further increased for 
women homozygous for the VDR haplotype 1 to 2.6 times the risk for fracture 
compared to women in ths reference group (Table 3a). When we analysed by type of 
fracture we observed the ADR genotype effect to be similar for vertebral fracture cases 
(p=0.07) and non-vertebral fracture cases (p=0.04;data not shown). The relative risk of 
fracture did not essentially change after adjustment for potential confounding factors 
such as age, weight, and bone density in the multivariate regression analysis. 

In this group of women we also determined the distribution of fractures according to 
COLIA1 genotype (Table 3b). In correspondence with what we previously found 10 
we observed the COLIA1 T allele to be associated with increased fracture risk, 
independent of BMD. Tb assess whether there was interaction between the VDR 
haplotype effect and thelCOLIAl genotype effect on fracture we determined the 
distribution of fractures According to VDR haplotype 1 in the different COLIA1 
genotype groups (Table 41. The distribution of fracture cases according to VDR 
genotype did not differ inuhe group of women with the COLIA1 GG genotype. 
However, in the COLIA1 risk groups of women with the GT and TT genotypes the 
distribution of fracture cases was strongly VDR genotype dependent (Table 4a). 
Logistic regression analysis snowed that the effect of VDR genotype on fracture risk is 
absent in women with the CQLIA1 GG genotype while the VDR genotype effect is 
large in the COLIA1 heterozygous GT and homozygous TT risk group (Table 4b). 
When age, VDR genotype, COLIAl genotype and fracture were considered together 
in a multivariate regression moHel we found that VDR genotype significantly modified 



# • 

WO 00/15839 PCT/EP99/07719 

15 

the COLIA1 gertotype effect (p=0.03 for the interaction term). The effect was found to 
be similar for noWertebral fracture cases and vertebral fracture cases and when bone 
density was enterad into the model the results did not change indicating the interaction 
effect to be independent of bone density. 

Although VDR gene polymorphisms have been implicated in the genetic regulation of 
BMD 7 , a meta-analysis showed that the effects on BMD are small 8 as we also 
demonstrated earlier in our study population 11 . While an early unpublished report 
j suggested VDR to predict osteoporotic fracture in older Australian women 17 , no. 
1 associations with osteoporotic fracture have been reported in two studies published 
/ until now on the relation of VDR with fracture 1 8 * 19 . These studies were small, 
I however, and used only the BsmI RFLP in their analysis. Our findings suggest the 
VDR to be involved in bone metabolic pathways other than those reflected in BMD 
but still leading to increased fracture risk. The VDR genotype dependent increased 
fracture risk is especially pronounced in interaction with COLIA1 genotype of which 
we already reported the COLIA1 Spl "T" allele to increase fracture risk in our study 
population 10 . The epistatic interaction between VDR and COLIA1 genotype we here 
describe points to biological interactions of the gene products. The VDR is a member 
of the steroid transcription factors, known to be important regulators of gene 
expression. Vitamin D dependent regulation of expression of bone-specific genes, 
such as osteocalcin, has been well documented 22 and also includes regulation of the 
expression of the collagen type I 1 at the level of transcription 20 * 21 . In RT-PCR 
experiments the Spl polymorphism has been shown to lead to differential binding 
affinity of the Spl transcription factor 9 and also to genotype dependent COLIA1 
mRNA expression levels 23 . Therefore, while the exact molecular mechanism 
underlying the associations we here describe remains to be elucidated, the VDR 
regulated expression of the collagen type I 1 gene could be an important factor in the 
synergistic interaction we observed and may differ across VDR alleles. 
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It is likely that interactions between genetic loci involved in a complex trait are a 
common phenomenon and several examples have already been demonstrated but 
mostly in model organisms. Our data represent the first example of interlocus 
interaction in relation to fractures between two well known candidate genes in 
osteoporosis, a complex trait in humans. We show that the interaction leads to 
increases in the risk of fracture, the clinically most relevant feature of osteoporosis, 
and that this increase in risk is independent of BMD, the most widely used diagnostic 
criterium for osteoporosis. This has important consequences not only for the analysis 
of the genetic basis of osteoporosis but also for the identification of individuals at risk 
of the disease. Finally, it also raises the possibility of developing new therapeutic 
intervention strategies based on the known involvement of the VDR and the COLIA1 
gene in bone metabolism. 

Methods 

Study subjects. The Rotterdam Study is a population-based cohort study of 7983 
subjects aged 55 or more years, residing in the Ommoord district of the city of 
Rotterdam in the Netherlands. The study was designed to document the occurrence of 
disease in the elderly in relation to several potential determinants 24 . A total of 10,275 
persons, of whom 9161 (89 percent) were living independently, were invited to 
participate in tht study in 1991. In the independently living subjects, the overall 
response rate was 77 percent for home interview and 71 percent for examination in a 
research centre, including measurement of anthropometric characteristics, bone 
mineral density and blood sampling. The Rotterdam Study was approved by the 
Medical Ethics Committee of the Erasmus University Medical School and written 
informed consent was obtained from each subject. The analysis of the association 
between COLIA1 genotype, VDR genotype and osteoporotic fracture was performed 
in a subgroup of women participating in the study. Baseline measurements of bone 
mineral density were available for 5931 independently living subjects from the study, 
but 1453 of these were excluded on the basis of age (>80 yrs), use of a walking aid, 
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diabetes mellitus or use of diuretic, estrogen, thyroid hormone or cytostatic drug 
therapy. From the 4478 remaining subjects, we studied a random sample of 1500 
women aged 55 to 80 years. Anthropometric data, DNA samples or genotype data for 
both loci were not available for 481 women, and we excluded women with the rare 
VDR haplotypes 4 and 5 (n=15) resulting in a final study group of 1004 women. 

Measurements. Height and weight were measured at the initial examination. Bone 
mineral density (in g/cxrfi) was determined by Dual Energy X-ray Absorptiometry 
(Lunar DPX-L densitometer; Lunar Corporation, Madison, WI, USA) at the femoral 
neck and lumbar spine (vertebrae L2 to L4) as described elsewhere. 25 Dietary intakes 
of calcium (mg/day) during the preceding year were assessed by food frequency 
questionnaire and adjusted for energy intake. Age at menopause and current cigarette 
smoking were assessed by questionnaire. For 732 women (73 percent), lateral 
radiographs of the spine from the fourth thoracic to the fifth lumbar vertebrae were 
obtained at baseline examination and analysed for the presence of prevalent vertebral 
fractures by morphometric analysis as previously described 26 The occurrence of 
incident non-vertebral fractures, including hip, wrist and other fractures, were 
recorded, confirmed and classified by a physician over a mean follow-up period of 3.8 
years. In total, 49 prevalent vertebral fracture cases and 52 incident non-vertebral 
fracture cases were recorded (7 hip, 6 upper humerus, 22 wrist, 4 hand, 4 ankle, 3 
foot, and 5 other fractures). Four subjects, in which both a vertebral and a nonvertebral 
fracture were present, were each counted as one fracture case, resulting in 97 cases 
with one or more fractures. 

Determination of COLIA1 and VDR genotypes. Genomic DNA was extracted from 
peripheral venous blood samples according to standard procedures and the 
polymorphism in the COLIA1 gene was detected by PGR with a mismatched primer 
that introduces a di-allelic restriction site, as previously described. 9 The test 
discriminates two alleles named S and s, corresponding to nucleotides G and T, 
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